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Abstract 

We study driven by an external electric field quantum orbital and spin dynamics of electron in a one-dimensional 
double quantum dot with spin-orbit coupling. Two types of external perturbation are considered: a periodic field 
at the Zeeman frequency and a single half-period pulse. Spin-orbit coupling leads to a nontrivial evolution in the 
spin and orbital channels and to a strongly spin- dependent probability density distribution. Both the interdot 
tunneling and the driven motion contribute into the spin evolution. These results can be important for the design 
of the spin manipulation schemes in semiconductor nanostructures. 
PACS numbers: 73.63.Kv,72.25.Dc,72.25.Pn 



Introduction 

Quantum dots, being one of the most intensively stu- 
died examples of natural and artificial nanostructures, 
attract attention due to the richness in the properties 
they demonstrate in the static and dynamic regimes [1]. 
A possible realization of qubits for quantum information 
processing can be done by using spins of electrons in 
semiconductor quantum dots [2], Spin- orbit coupling 
makes the dynamics even in the basic systems such as 
the single-electron quantum dots extremely rich both in 
the orbital and spin channels. If the frequency of the 
electric field driving the orbital motion matches the 
Zeeman resonance for electron spin in a magnetic field, 
the spin-orbit coupling causes a spin flip. This effect 
was proposed in refs. [3,4] to manipulate the spin states 
by electric means. The efficiency of this process is much 
greater than that of the conventional application of a 
periodic resonant magnetic field. The ability to cause 
coherently the spin flip in GaAs quantum dots was 
demonstrated in ref. [5] where the gate-produced elec- 
tric field induced the spin Rabi oscillations. In ref. [6] 
periodic electric field caused the spin dynamics by indu- 
cing electron oscillations in a coordinate-dependent 
magnetic field. In addition, these results confirmed that 
the spin dephasing in GaAs quantum dots, arising due 
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to the spin-orbit coupling [7,8] is not sufficiently severe 
to prohibit a coherent spin manipulation. 

The spin dynamics experiments [5,6] necessarily use at 
least a double quantum dot to detect the driven spin 
state relative to the spin of the reference electron. Multi- 
ple quantum dots realizations become nowadays the sub- 
ject of extensive investigation [9]. In double quantum 
dots an interesting charge dynamics occurs and requires 
theoretical understanding. In this article we address full 
driven by an external electric field spin and charge quan- 
tum dynamics in a one-dimensional double quantum dot 
[10-13]. Despite the simplicity, these systems show a rich 
physics. In the wide quantum dots, where the tunneling 
is suppressed, and the motion is classical, the interdot 
transfer occurs only due to the over-the-barrier motion, 
and a chaos-like behavior is usually expected. The irregu- 
lar driven behavior in the spin and charge dynamics in 
these systems was studied in ref. [14]. In the quantum 
double quantum dots, the tunneling between single 
quantum dots is crucial and the spin-orbit coupling 
makes the interdot tunneling spin-dependent [15-17]. In 
quantum systems a finite set of energy eigenstates allows 
only for a strongly irregular rather than a real chaotic 
behavior. These orbital and spin dynamical irregularities 
are important for the understanding of the quantum pro- 
cesses in multiple quantum dots. 

In this article we consider various regimes for a one- 
dimensional double quantum dot with spin-orbit cou- 
pling driven by an external electric field and analyze the 
probability and spin density dynamics in these systems. 



SpringerOpen 0 



© 2011 Khomitsky and Sherman; licensee Springer. This is an Open Access article distributed under the terms of the Creative 
Commons Attribution License (http://creativecommons.0rg/licenses/by/2.O), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Khomitsky and Sherman Nanoscale Research Letters 2011, 6:212 
http://www.nanoscalereslett.eom/content/6/1/212 



Page 2 of 5 



Hamiltonian, time evolution, and observables 

We use a quartic potential model to describe a one- 
dimensional double quantum dot [18], 



17(f) = U 0 {-2(x/d) 2 + {x/df), 



(1) 



where the minima located at d and -d are separated 
by a barrier of height U 0 , as shown in Figure 1. We 
assume that the interminima tunneling is sufficiently 
weak such that the ground state can be described with a 
high accuracy as even linear combination of the oscilla- 
tor states with a certain "harmonic" frequency a>o 
located near the minima. The double quantum dot is 
located in a static magnetic field B z along the z-axis and 
is driven by an external electric field S'(t) parallel to the 
#-axis. The full Hamiltonian H = Hq + H so + V , where 
the time-independent parts are given by 



Hn 



2m 



U{x) 



A, 



H so = [Pcr x + ao y )p x , 
and the time-dependent perturbation is 

V = -e£{t)x. 



(2) 



(3) 



(4) 



Here p x is the momentum operator, m is the electron 
effective mass, e is the electron charge, A z = \g\fisB z (we 
assume below g < 0) is the Zeeman splitting, and <7; are 
the Pauli matrices. The electron Lande factor g deter- 
mines the effect of B z , which in this geometry is reduced 
to the Zeeman spin splitting only. The bulk-originated 
Dresselhaus (J3) and structure-related Rashba (a) para- 
meters determine the strength of spin-orbit coupling 
and make the electron velocity defined as 



v = x= t[ h o + H so , x] =p x /m + fio x + aa y , 



(5) 



spin-dependent. 
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Figure 1 A schematic plot of the double-well potential 
described by Equation (1). Double green (red) lines correspond to 
the spin-split even (odd) tunneling-determined orbital states. 



We use the highly numerically accurate approach to 
describe the dynamics with the sum of Hamiltonians in 
Equations (2)-(4). As the first step we diagonalize 
exactly the time-independent H 0 + H so in the truncated 
spinor basis y/ n (x) |<t) of the eigenstates of the quartic 
potential in magnetic field without spin-orbit coupling 
with corresponding eigenvalues E na . As a result, we 
obtain the basis set |(// n ) where bold n incorporates the 
spin index. For the presentation, it is convenient to 
introduce the four-states subset: \y/{) = y/i(#)|T), ^2) = 
Vi(x)\l), \y/ 3 ) = y/2M|f)> \W*> = V2WII). and to note 
that the spin-dependent bold index may not correspond 
to the state energy due to the Zeeman term in the 
Hamiltonian. The wavefunction y/\{x) ((//2M) is even 
(odd) with respect to the inversion of x. In the case of 
weak tunneling, assumed here, these functions can be 
presented in the form: ^1,2 W = (VtM ± ^W)/V2, 
where y/^{x) and y/B.(x) are localized in the left and in 
the right dot, respectively. 

As the second step we build in the full basis the 
matrix of time-dependent V and study the full 
dynamics with the wavefunctions: 



(6) 



The expansion coefficients C„(f) are then calculated as: 
d 



dt 



§n(t) = ij-£(t) ^2%m{t)x nm e~ 



i{E m -E n )t/h 



(7) 



Where *nm = (fn\x\ f m ) . The spin-dependence of the 
matrix element of coordinate responsible for the spin 
dynamics is determined with 



i(E n - E m )x nm = h(i/ n \i>\^ n 



(8) 



and the spin-dependent velocity in Equation (5). 

With the knowledge of the time-dependent wavefunc- 
tions (6) one can calculate the evolution of probability 
p(x, t) and spin Si(x, f)-density 



p{x,t) = W(x,t)9{x,t), 



Si[x, t) = * t (x, t)o-;*(x, t). 



(9) 



(10) 



Since we are interested in the interdot transitions, 
with these distributions we find the gross quantities, e. 
g., for the right quantum dot: 



Jo 

/>oo 

°r(0 = / Si(x, t)dx, 
Jo 



(11) 



(12) 



Khomitsky and Sherman Nanoscale Research Letters 2011, 6:212 
http://www.nanoscalereslett.eom/content/6/1/212 



Page 3 of 5 



where a) R (t) is the probability to find electron and 
is the analog of expectation value of the spin 
component. 

Calculations and results 

As the electron wavefunction at t = 0 we take linear 
combinations of two out of four low-energy states. The 

initial state in the form {ir^x) ± f 2 {x))\ \) j \/I is 

localized in the left quantum dot, corresponding to the 

parameters = fs(0)l / 

Two types of electric field were considered as the 
external perturbation. The first one is the exactly peri- 
odic perturbation for all t > 0: 

£{t) = £ 0 sin (27rt/T z (B z )), (13) 

Where T Z (B Z ) = 2nh/A z is the Zeeman period. The 
second type is a half-period pulse, same as in Equation 
(13), but acting at the time interval 0 <t <T Z {B Z )I2 only. 
The spectral width of the pulse covers both the spin 
and the tunneling splitting of the ground state, thus, 
driving the spin and orbital dynamics simultaneously. 
Since Tz(Bz)w » 1, that is the corresponding frequen- 
cies are much less than those for the transitions 
between the orbital levels corresponding to a single dot, 
the higher-energy states follow the perturbation adiaba- 
tically. The field strength S 0 is characterized by para- 
meter / such that \e\<% 0 = f x U 0 /2d. Here we 
concentrate on the regime of a relatively weak coupling 
(/=» 1). 

Where the shape of the quartic potential remains 
almost intact in time, and the interdot tunneling is still 
crucially important. For the magnetic field we consider 
two different regimes A z = AE g /2 and A z = 2AE g to illus- 
trate the role of the Zeeman field for the entire 
dynamics. 

We consider a nanostructure with d = 25^/2 nm and 
U Q = 10 meV. The four lowest spin-degenerate energy 
levels are E 1 = 3.938 meV, E 2 = 4.030 meV, E 3 = 
9.782 meV, £ 4 = 11.590 meV counted from the bottom 
of a single quantum dot with the tunneling splitting AE g 
= E 2 - Ei = 0.092meV, and the corresponding timescale 
2nh/AE g = 45ps. The spin-orbit coupling is described by 
parameters a = 1.0 • 10~ 9 eVcm and /3 = 0.3 • 10~ 9 
eVcm. The field parameter /= 0.125, corresponding to 
<S?o = 177 V/cm. We use the truncated basis of 20 states 
with the energies up to 42 meV. 

We begin with the exactly periodic driving force, as 
illustrated in Figure 2 where |f n | 2 for three states are 
presented. Since the motion is periodic, here we use the 
Floquet method [13,19,20] based on the exact calcula- 
tion at the first period and then transformed into the 
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Figure 2 Motion driven by the exactly periodic field. Upper 
panel: B z = 1.73T, A z = AF 9 /2 and T z (fij = 90 ps; lower panel: B z = 
6.92T, A z = 2A£g, and T Z (B Z ) = 22 ps. The states for c n (f) are marked 
near the plots. The upper panel demonstrates a relatively slow 
dynamics on the top of the fast oscillations. The increase in the £ 2 (f) 
term corresponds to the possible spin-flip due to the external 
electric field. 



integer number of periods. Figure 2 demonstrates the 
interplay between the tunneling and the spin-flip pro- 
cess. The results indicate that the exact matching of the 
driving frequency with the Zeeman splitting generates 
the spin flip which is clearly visible as the initial spin-up 
(fi and f 3 ) components are decreasing to zero and, at 
the same time, the opposite spin-down components (f 2 
and Cn) reach their maxima (not shown in the upper 
panel). The spin-flip time is approximately 350 T z (B z ) 
(or 31 ns) for the weak magnetic field (upper panel) and 
24T Z (B z ) (or 528 ps) for the strong field (lower panel). 
Such an increase in the Rabi frequency with increasing 
magnetic field is consistent with previous theoretical 
[3,4] and experimental results [5]. 
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Figure 3 Motion driven by the single-pulse field. Upper panel: 
B z = 1.73T, A z = Afg/2, and 7 Z (B Z ) = 90 ps; lower panel: B z = 6.92T, 
A z = 2AF g , and f z (6 z ) = 22 ps Black line is the probability to find the 
electron in the right quantum dot. Red and blue dashed lines show 
corresponding spin components, as marked near the lines, 
determined both by the spin-orbit coupling and external magnetic 
field. 



of x -operator for the pairs of states with opposite space 
parity having the same dominating spin projection. 
Hence, the dynamics involving all four lowest levels first 
of all triggers the transitions inside these pairs which do 
not involve the spin flip and only after this the spin-flip 
processes can become significant. As a result, Figure 3 
shows that the spin flip has only partial character while 
the free tunneling dominates as soon as the pulse is 
switched off. A detailed description of other processes 
of nonresonant driven dynamics in the case of a half- 
period perturbation can be found in ref. [21]. 

Conclusions 

We have studied the full driven quantum spin and 
charge dynamics of single electron confined in one- 
dimensional double quantum dot with spin-orbit cou- 
pling. Equations of motion have been solved in a finite 
basis set numerically exactly for a pulsed field and by 
the Floquet technique for the periodic fields. We 
explored here the regime of relatively weak coupling to 
the external field, where a nontrivial dynamics already 
occurs. Our results are important for the understanding 
of the effects of spin-orbit coupling for nanostructures 
as we have demonstrated a possibility to achieve a con- 
trollable spin flip at various time scales and in various 
regimes by the electrical means only. 
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As the second example we consider the probabilities 
ta R (i) and Ofl(t) for the pulse-driven motion, presented 
in Figure 3. As one can see in the figure, the initial 
stage is the preparation for the tunneling, which devel- 
ops only after the pulse is finished. Electric field of the 
pulse induces the higher-frequency motion by involving 
higher-energy states, as can be seen in the oscillations at 
t < T z (B z )/2, however, prohibits the tunneling. Such a 
behavior of the probability and spin density can be 
explained by taking into account the detailed structure 
of matrix elements x nm . Namely, due to the symmetry 
of the eigenfunctions in a symmetric double QW the 
largest amplitude can be found for the matrix element 
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